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We introduce a coarse-grained description of the surface energy of a nematic liquid crystal. The 
thermal fluctuations of the nematic director close to the surface renormalize at macroscopic scales the 
bare surface potential in a temperature-dependent way. The angular dependence of the renormalized 
potential is dramatically smoothed, thus explaining the success of the Rapini-Papoular form. Close 
to the isotropic phase, the anchoring energy is strongly suppressed and the change of its shape allows 
for anchoring transitions. Our theory describes quantitatively the temperature dependence of the 
anchoring energy and the temperature-induced anchoring transitions reported in the literature. 
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In recent years, surface phenomena have attracted 
a lot of interest. Particularly, the interface between 
liquid crystals jjj and solid substrates displays a rich 
variety of behaviors: orientational wetting and 
spreading HQ, memory effects M, surface melting ||, 
Kosterlitz-Thouless transitions lid], quasi-critical be- 
havior of surface energies Jll|-|l3|], surface anchoring 
transitions 14-22], etc. Despite the complexity and 
the diversity of the interactions between nematic liq- 
uid crystals and solid substrates, some simple — though 
unexplained — ubiquitous aspects emerge from the exper- 
iments: the angular dependences of the surface potentials 
are extremely smooth and well described by the so-called 
Rapini-Papoular law |23| , at high temperatures the pre- 
ferred surface orientation is usually either parallel or per- 
pendicular to the substrate, tilted orientations are diffi- 
cult to achieve p4| , and temperature-driven anchoring 
transitions systematically occur immediately below the 
bulk isotropic transition temperature |2dJ. In this Let- 
ter, we show that all these effects can be explained in 
terms of a renormalization of the surface energy by the 
short-wavelength orientational fluctuations in the bulk. 

Nematic liquid crystals are fluid mesophases made of 
elongated molecules displaying a broken orientational 
symmetry along a non-polar direction called the nematic 
director n |Q. At a molecular level, nematics exhibit 
large orientational fluctuations and usually some degree 
of short-range biaxial and positional order. The link be- 
tween the microscopic and the macroscopic description 
can be established by means of a coarse-graining proce- 
dure |25|] . Before determining the consequences of this 
procedure on the surface behavior, let us briefly discuss 
how it is carried out in the bulk. One can associate a local 
orientation m to each molecule; then the probability p[m] 
of a given instantaneous microscopic configuration m(r) 
is proportional to the Boltzmann factor exp(— /3.F[m]), 
where ^"[m] is a complex microscopic free-energy ]2(| 
and (3 = is the inverse temperature. At macro- 

scopic scales the bulk free-energy F[n] is well described 



by a simple elasticity involving only the average molec- 
ular orientation n and its gradient [^,0. The quantity 
exp(— /Si^n]) gives the probability of observing a given 
smooth director configuration n(r), whatever the details 
of its rapidly varying components. One therefore obtains 
the coarse-grained free-energy F[n] by summing the mi- 
croscopic probabilities exp(— /3.F[m]) over all the Fourier 
components m(q) with wavevector |q| > A, where A 
is some macroscopic cutoff, the components m(q) with 
|q| < A being fixed and denoted by n(q). The resulting 
free-energy F[n], which corresponds to the Frank elas- 
ticity p7| ], is therefore meaningful only for the slowly 
varying Fourier components |q| < A of n. Such a coarse- 
graining procedure yields a renormalized elasticity ex- 
panded in power series of the derivatives of n, whose 
coefficients depend on T and are expected to scale as 
powers of the range b of the molecular interactions. This 
is why, for Kb <C 1, it is justified to retain in F[n] only 
the lowest-order gradient terms. 

Using the Frank elasticity F[n] implicitly entails a 
coarse-graining in the bulk. For the sake of consistency, 
this procedure should also be performed on the surface. 
The latter will then be effectively transformed into a 
blurred layer of width A -1 acquiring some properties of 
the bulk. Such a coupling between surface and bulk is 
usually introduced in a phenomenological way by means 
of Landau expansions in the tensorial nematic order- 
parameter Q P |^J^- [T^ , p8| , p9t . However, modeling sur- 
face properties in this way is rather complicated, since 
one has to deal with spatially varying tensorial fields, 
and also somewhat arbitrary, since high powers of Q 
should be included due to the first-order character of the 
nematic-isotropic transition. 

We have calculated the renormalization of the surface 
energy by coarse-graining the director orientation from a 
mesoscopic length a, of the order of the nematic coher- 
ence length £ni at the isotropic transition |l| , up to a scale 
A" 1 at least a few times larger than a. Our results, based 
on perturbation theory, are valid for weak anchorings, in 
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FIG. 1. Schematic representation of the surface and its an- 
choring potential at the mesoscopic scale a and coarse-grained 
scale A" 1 . The mesoscopic potential v{9) and its correspond- 
ing effective coarse-grained potential V(9) are plotted in po- 
lar coordinates in arbitrary units. The renormalized potential 
V{6), which depends on temperature, favors here an oblique 
anchoring at Ti and an homeotropic one at T2 > Ti. 



the sense £ 3> £ni, where £ is the anchoring extrapolation 
length 0. Calling and 4> the director's polar angles 
and starting from a "bare" surface anchoring energy ex- 
panded in Fourier harmonics of the form cs(2n9) cs(m(/>) 
(cs being either cos or sin), we find that each harmonic 
is independently renormalized by a Debye- Waller factor 
exp(— a n — (3 m ), with a n ,(3 n oc n 2 k^T/Ka, where K is a 
bulk nematic elastic constant. At scales larger than A -1 , 
the high-order harmonics are thus strongly suppressed 
by the director thermal fluctuations: this explains the 
success of the Rapini-Papoular form |22(| . Moreover, the 
anchoring energy naturally acquires a temperature de- 
pendence through the elastic constant K{T). The sur- 
face energy is thus reduced close to the nematic-isotropic 
transition — where K is lowered — in agreement with ex- 
periments ]TT| [l3[] . The different temperature depen- 
dences of the surface harmonics allow for anchoring tran- 
sitions: as the temperature increases, the suppression of 
the high-order harmonics shifts the minimum of the an- 
choring energy towards some symmetry axis of the sur- 
face. Our results fit well the quasi-critical temperature 
dependence of the azimuthal anchoring energy measured 
by Faetti et al. pT| , the oblique-to-homeotropic anchor- 
ing transition measured by Patel and Yokoyama pfj| , and 
the bistable oblique-to-planar anchoring transition fit] ] 
measured by Jagemalm et al. |2^| . It turns out that the 
effective macroscopic anchoring, whose minimum gives 
the coarse-grained director orientation at the surface, can 
significantly differ from the mesoscopic surface potential 
(Fig. |l|). The director fluctuations within the coarse- 
grained region dramatically smooth the fine details of the 
anchoring energy. They can also shift the average equi- 
librium position of the director at the surface, similarly 
to the amplitude-dependent shift of the average position 
of an anharmonic oscillator [ROf . 



Precisely, we consider a semi-infinite nematic slab in 
the z > half-space and we describe the nematic director 
by its spherical coordinates 9, (f> centered on the z-axis. 
We start with a bulk elasticity already coarse-grained 
on a mesoscopic length a ss £ni, such that it can be 
expressed in the usual Frank form. In the one-constant 
approximation [jjj, its harmonic part, expanded about an 
arbitrary direction (0q, 4>o)i takes the form 



Ho = \K{T) J d 3 r [(V0) 



sm 2 9 (y^y 



(1) 



At the length-scale a, the "bare" surface potential H s is 
given by some local functional of the mesoscopic direc- 
tor's orientation at the surface 



H s = d 2 r ± v( 



(2) 



The total free-energy Ft of the nematic slab is given by 
the path integral 



exp(-/3F t ; 



V[9]V[4>] exp [-/3(H + W„ h + H.)\ 



(3) 



where 7i n h contains the non-harmonic bulk terms. To 
further coarse-grain on a length-scale A" 1 > a, we put 
0(r) = 0<(r) + 0>(r) and cf>(r) = 0<(r) + </>>(r), where 
< (r), < (r) have Fourier components with wavevectors 
|q| < A, and > (r), > (r) have Fourier components with 
A < |q| < 2ir/a. Integrating out 9 > and cj) > yields the 
renormalized Hamiltonian 7i < at length-scales A -1 ; to 
lowest-order in perturbation theory, it is given by |25| 



H < = Ha + (Wnh + H s 



(4) 



where (••■)> indicates the statistical average over the 
high-wavevector components, weighted by the Gaussian 
Hamiltonian Hq. At first-order, the bulk and the surface 
are therefore renormalized independently. 

Due to the n — > — n invariance, the bare surface energy 
density can be Fourier expanded as 



cos n0 [w^ m cos mcf) + w^ m sin mq 



- sin n9 [w s ^ m cos mcf> + w s ^ m sin m<j>] , (5) 



where n runs over all even integers, and m runs over all 
even (resp. odd) integers in the first (resp. second) sum. 

Even though is non- linear, the renormalized Hamil- 
tonian ^ can be transformed into Gaussian integrals by 
writing the trigonometric functions as complex exponen- 
tials. We obtain a renormalized surface energy density 
V(9 < ,(f) < ) whose Fourier components W"m (a, being 
cither c or s) are separately renormalized: 
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FIG. 2. Azimuthal ancho ring energy Wo vs. temperature 
T according to Faetti et al. jll]] . Points: experimental data; 
solid line: our fit. 



W«£(t)=<&exp 



where, for Aa <C 1, 



t = 



'W 2 )>^ B r/ 
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d 3 <? 1 k B T 



(2n) 3 Kq 2 ixKa 



(6) 



(7) 



Note that (0 > (r) 2 ) > = i/sin 2 9q also appears in (^). 

It follows from (g) that for large t, i.e., in the vicinity of 
the ncmatic-isotropic transition, where K is reduced, the 
effective surface potential should be very-well described 
by its first harmonics. This is indeed what was measured 
by Faetti et al. using the nematic liquid crystal 5CB, 
planarly anchored on a glass plate treated by oblique 
evaporation of SiO: exploring the azimuthal anchoring 
energy in regions far away from the parabolic minimum, 
they could not detect any deviation from the simple form 



V = W sin 2 tf> = W§2 cos 20 + const 



(8) 



Setting 9 = 
with flm, i.e. 



tt/2, </> = 0, we have fitted their data 



>exp(— 2t), in which we have 



assumed a Landau-de Gennes form for the temperature 
dependence of K j3l| , yielding 
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TrK m a \3 + V9 - 8? 



(9) 



where r = (T - T*)/(T OT - T*) is the reduced tem- 
perature, Tni the nematic-isotropic transition temper- 
ature, T* the isotropic supercooling temperature, and 
i^Ni the elastic constant at the transition. With Tni — 
306.8 °K, T N i - T* ~ 1.5 °K @, our best fit, shown 
in Fig. ||, yields *ni = k^T^i/irKm a = 2.1 ± 0.1 and 
W§%(t m ) = (-2.08±0.08)xl0~ 4 erg/cm 2 . For an extrap- 
olated Km ~ 0.5 pN [R3|, the value of tm yields a ~ 13 A, 
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FIG. 3. Surface easy-axis 8 (measured with respect to the 
surface normal) vs. temperature T according to Patel and 
Yokoyama Eu|. Points: experimental data; solid line: our 
fit. Inset: anchoring phase diagram (in arbitrary units) cor- 
responding to Eq. (LOI). 



which roughly compares with £ni for a first-order phase 
transition. 

Close to the nematic-isotropic transition of the ne- 
matic compound E7, Patel and Yokoyama [^0| have ob- 
served a tilted-to-homeotropic anchoring transition on a 
fluoropolymer-coated surface. Again, we can explain it 
by retaining in V the first two harmonics allowed by sym- 
metry, i.e., 



V = W^ cos 29 + W^ cos 49 



40 ' 



(10) 

Minimizing ( jToj ) with respect to 9 yields the anchoring 
phase diagram shown in the inset of Fig. ||, which displays 
the following regions separated by second-order transi- 
tion lines: conical oblique (O), homeotropic (H), degen- 
erated planar (P), and metastable homeotropic/planar 
(H/P). According to (|) and (§), as temperature in- 
creases, the easy-axis of the anchoring follows the typical 
path shown in the phase diagram. Setting Tni — 57.7 °C 
and Tjj = 54.48 °C §0), where T H is the homeotropic- 
to-oblique transition temperature, we fit closely the tem- 
perature dependence of the easy-axis direction 9 with 
T N i - T* ~ 1.62 °K and t m = 1.2 ± 0.04 (Fig. §). 

Using the same values of Tni, Tni — T* and iNi, which 
depend only on the nematic material, we can fit equally 
well the bistable oblique-to-planar anchoring transition 
observed by Jagemalm et al. Q at T = T P = 49.7 °C 
for the same compound E7 (FigTH) . This anchoring tran- 
sition also appears close to the nematic-isotropic tran- 
sition, in the presence of a glass substrate treated by 
an oblique evaporation of SiO. Here, by symmetry, the 
lowest-order expansion of V is 

V = W^a cos 26> + W C | cos 2<j> + WfJ sin 29 cos 0. (11) 

Setting We = W^/W£ and W B< p = \W%/W&\, the 
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FIG. 4. Surface easy-axis (i 
ing to Jagemalm et al. |^2|. Points: experimental data; solid 
line: our fit. Inset: anchoring phase diagram corresponding to 
Eq. (FLIT) ; the arrow indicates the path followed as T increases. 



corresponding anchoring phase diagram, shown in the 
inset of Fig. ||, displays the four regions, separated by 
second-order transition lines, that are observed experi- 
mentally p6|JT7| ] : oblique in the evaporation plane (O), 
bistable symmetric with respect to the evaporation plane 
(B), homeotropic (H), and planar orthogonal to the evap- 
oration plane (P). Once fixed the values of Tp, and, by 
the previous fit, of Tni, ?ni — T* , and tm, the fit on 9 
has no adjustable parameter, while for the fit on <fi there 
remains the only free parameter ifli/^oii best adjusted 
to the value 6.76 ± 0.05. 

To estimate the validity of our first-order perturbative 
expansion in the surface potential, we have calculated the 



second-order correction to W%q coming from w 
Aa <§; 1, we find 
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which is negligible for A -1 <C I, where I = K/w^q is a 
bare anchoring extrapolation length; this sets the limit 
of validity of our analysis. Coarse-graining on a length 
A -1 > £ effectively reduces the anchoring for purely elas- 
tic reasons, as will be described elsewhere. 

Finally, note that, since we coarse-grained on a length 
A -1 > £ni, our model is insensitive to the variations of 
the scalar order-parameter S This does not imply, 
however, that S must be uniform throughout the sam- 
ple: any substrate inducing a surface variation of S will 
still be described — at length-scales larger than £ NI — by 
some potential v{6, <p) depending only on the director's 
orientation. Our analysis simply disregards the underly- 
ing surface variations of S: this is a limitation only when 
£ni becomes critically large. 
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